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Abstract: At the southern site of the Pierre Auger Observatory, which is close to completion, an exposure 
that significantly exceeds the largest forerunner experiments has already been accumulated. We report 
a measurement of the cosmic ray energy spectrum based on the high statistics collected by the surface 
detector. The methods developed to determine the spectrum from reconstructed observables are described. 
The energy calibration of the observables, which exploits the correlation of surface detector data with 
fluorescence measurements in hybrid events, is presented in detail. The methods are simple and robust, 
exploiting the combination of fluorescence detector (FD) and surface detector (SD) and do not rely on 
detailed numerical simulation or any assumption about the chemical composition. Besides presenting 
statistical uncertainties, we address the impact of systematic uncertainties. 
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The Pierre Auger Observatory [1] is designed to 
measure the extensive air showers produced by the 
highest energy cosmic rays (E > 10 18 5 eV) with 
the goal of discovering their origins and composi- 
tion. Two different techniques are used to detect air 
showers. Firstly, a collection of telescopes is used 
to sense the fluorescence light produced by exci- 
tation of nitrogen induced by the cascade of parti- 
cles in the atmosphere. The FD provides a nearly 
calorimetric, model-independent energy measure- 
ment, because the fluorescence light is produced 
in proportion to energy dissipation by a shower in 
the atmosphere [2, 3]. This method can be used 
only when the sky is moonless and dark, and thus 
has roughly a 10% duty cycle [4]. The second 
method uses an array of detectors on the ground to 
sample particle densities as the air shower arrives 
at the Earth's surface. The surface detector has a 
100% duty cycle [5]. A subsample of air show- 
ers detected by both instruments, dubbed hybrid 
events, are very precisely measured [6] and pro- 
vide an invaluable energy calibration tool. Hybrid 
events make it possible to relate the shower energy 
(FD) to the ground parameter 5(1000). 



The parameter 5(1000) characterises the energy of 
a cosmic ray shower detected by the SD array and 
is the signal in units of VEM that would be pro- 
duced in a tank at a distance of 1000 m from the 
shower axis. One VEM is the signal produced by a 
single relativistic muon passing vertically through 
the centre of a water tank. A likelihood method is 
applied to obtain the lateral distribution function, 
where the shower axis, 5(1000) and the curvature 
of the shower front are determined [7]. The se- 
lection criteria are such to ensure the rejection of 
accidental triggers (physics trigger) and the events 
are well contained in the SD array (quality trig- 
ger), i.e. we require that all six nearest neighbours 
of the station with the highest signal be active. In 
this way we guarantee that the core of the shower 
is contained inside the array and enough of the 
shower is sampled to make an S(1000) measure- 
ment. The present data set is taken from 1 January, 
2004 through 28 February, 2007 while the array 
has been growing in size. To ensure an excellent 
data quality we remove periods with problems due 
to failures in data acquisition, due to lightning and 
hardware difficulties. We select events only if the 
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Figure 1 : Integral number of events vs cos 2 9 for 
the indicated minimum value of 5(1000). 

zenith angle is less than 60° and the reconstructed 
energy is above 3 EeV. For this analysis, the array 
is fully efficient for detecting such showers, so the 
acceptance at any time is solely determined by the 
geometric aperture of the array [8]. The integrated 
exposure mounts up to about 5165 km 2 sr yr, 
which is a factor of more than 3 larger than the 
exposure obtained by the largest forerunner ex- 
periment AGASA [9]. Moreover the present ac- 
ceptance exceeds the one given in [10] by a fac- 
tor of about 3. For a given energy the value of 
5(1000) decreases with zenith angle, 6, due to at- 
tenuation of the shower particles and geometrical 
effects. Assuming an isotropic flux for the whole 
energy range considered, i.e. the intensity distribu- 
tion is uniform when binned in cos 2 9, we extract 
the shape of the attenuation curve from the data. In 
Figure 1 several intensities, I{ = I(> 5;(1000)), 
above a given value of lg 5^(1000) are shown as 
a function of cos 2 9. The choice of the threshold 
lg 5(1000) is not critical since the shape is nearly 
the same within the statistical limit. The fitted at- 
tenuation curve, CIC(9) = 1 + a x + b x 2 , is a 
quadratic function of x — cos 2 9 — cos 2 38° as dis- 
played in Figure 2 for a particular constant inten- 
sity cut, To = 128 events, with a = 0.94 ± 0.06 
and b = — 1.21 ± 0.27. The cut corresponds 
to a shower size of about 538° = 47VEM and 
equivalently to an energy of about 9 EeV. Since 
the average angle is (9) ~ 38° we take this an- 
gle as reference and convert 5(1000) into 538° by 
5 38 o = S (1000) /CIC {9). It may be regarded 
as the signal 5(1000) the shower would have pro- 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cos 2 e 

Figure 2: Derived attenuation curve, CIC(9), fit- 
ted with a quadratic function. 

duced had it arrived at 9 — 38°. The reconstruc- 
tion accuracy of the parameter 5(1000), 05(1000)' 
comprises 3 contributions and these are taken into 
account in inferring 538° an d its uncertainty <Js ago '■ 
a statistical uncertainty due to the finite size of the 
detector and the limited dynamic range of the sig- 
nal detection, a systematic uncertainty due to the 
assumptions of the shape of the lateral distribu- 
tion and finally due to the shower-to-shower fluc- 
tuations [11]. To infer the energy we have to es- 
tablish the relation between 538° and the calori- 
metric energy measurement, Efd- A set of hybrid 
events of high quality is selected based on the crite- 
ria reported in [6] without applying the cut on the 
field of view, which appears to have a negligible 
effect for the topic addressed here. A small correc- 
tion to account for the energy carried away by high 
energy muons and neutrinos, the so-called invisi- 
ble energy, depends slightly on mass and hadronic 
model. The applied correction is based on the av- 
erage for proton and iron showers simulated with 
the QGSJet model and sums up to about 10% and 
its systematic uncertainty contributes 4% to the to- 
tal uncertainty in FD energy [3]. Moreover the SD 
quality cuts described above are applied. The cri- 
teria include a measurement of the vertical aerosol 
optical depth profile (VAOD(h)) [12] using laser 
shots generated by the central laser facility (CLF) 
[13] and observed by the FD in the same hour of 
each selected hybrid event. The selected hybrid 
events were used to calibrate the SD energy. The 
following procedure was adopted. For each hy- 
brid event, with measured FD energy Efd, the 
SD energy estimator 538° was determined from the 
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Figure 3: Correlation between lg Efd and lg S^go 
for the 387 hybrid events used in the fit. The full 
line is the best fit to the data. Events below the 
dashed line were not included in the fit. 



measured ,5(1000) by using the constant intensity 
method described above. For each event the uncer- 
tainty in S380 is estimated by summing in quadra- 
ture three contributions: the uncertainty in the con- 
stant intensity parametrization, <Js 3S o (CJC) , the 
angular accuracy of the event, a cos g, and the uncer- 
tainty in the measured S'(IOOO), crg(iooo)- The flu- 
orescence yield used to estimate the energy Efd is 
taken from [14]. An uncertainty in the FD energy, 
&e fd , was also assigned to each event. Several 
sources were considered. The uncertainty in the 
hybrid shower geometry, the statistical uncertainty 
in the Gaisser-Hillas fit to the profile of the en- 
ergy deposits and the statistical uncertainty in the 
invisible energy correction were fully propageted. 
The uncertainty in the VAOD measurement was 
also propagated to the FD energy on an event- 
by-event basis, by evaluating the FD energy shift 
obtained when changing the VAOD profile by its 
uncertainty. These individual contributions were 
considered to be uncorrected, and were thus com- 
bined in quadrature to obtain o\e fd . The data 
appear to be well described by a linear relation 
IgEpo = A + B ■ lg (see Figure 3). A lin- 
ear least square fit of the data was performed. To 
avoid possible biases, low energy events, below the 
dashed line, which is orthogonal to the best fit line 
and intersects it at lg(5 , 38<= = 15 VEM), were not 
included in the fit. 
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Figure 4: Fractional difference between the FD and 
SD energy for the 387 selected hybrid events. 

An iterative procedure was used to determine the 
dashed line, and it was checked that the results 
of the fit were stable. The best fit yields A = 
17.08 ± 0.03 and B = 1.13 ± 0.02 with a reduced 
X 2 of 1.3 for lg £5,0 = A+B-lgS 3 &° in [eVJ.The 
relative statistical uncertainty in the derived SD en- 
ergy, <j e sd /e S d, is rather small, e.g. of the order 
of 5% at 10 20 eV. The energy spectrum J is dis- 
played in Figure 5 together with its statistical un- 
certainty. The individual systematic uncertainties 
in determining Esd coming from the FD sum up 
to 22%. For illustrative purposes we show in Fig- 
ure 6 the difference of the flux with respect to an 
assumed flux oc E~ 2M . The largest uncertainties 
are given by the absolute fluorescence yield (14%), 
the absolute calibration of the FD (9.5%) and the 
reconstruction method (10%). The uncertainty due 
to the dependence of the fluorescence spectrum 
on pressure (1%), humidity (5%) and temperature 
(5%) are take into account as well as the wave- 
length dependent response of the FD, the aerosol 
phase function, invisible energy and others, which 
are well below 4% (see [4] for details). 

Discussion and outlook 

When inferring the energy spectrum from SD data 
we utilise the constant intensity method to cali- 
brate the SD data. The systematic uncertainties 
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Figure 5: Auger spectrum J 
as a function of energy. Ver- 
tical error bars represent the 
statistical uncertainty only. 
The statistical and system- 
atic uncertainties in the en- 
ergy scale are of the order of 
pa 6% and ps 22%, respec- 
tively. 
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have been scrutinised and the resulting spectrum 
is given. Several activities are on-going to reduce 
the systematic uncertainties of the energy estimate, 
e.g. the detector calibration uncertainty and the un- 
certainty of the fluorescence yield. Reducing these 
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Figure 6: Fractional difference between the de- 
rived spectrum and an assumed flux oc E~ 2M as 
a function of energy. 

uncertainties will make it desirable to deconvolve 
the energy spectrum using the estimate of the en- 
ergy resolution. The presented spectrum is com- 
pared with a spectrum derived on basis of hybrid 
data only in T. Yamamoto et al. [15]. Astrophysi- 
cal implications are also discussed there. 
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